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Emission of bremsstrahlung photons accompanying proton decay of nuclei is studied. The new improved 
multipolar model describing such a process is presented. The angular formalism of calculations of the matrix 
elements is stated in details. The bremsstrahlung probabilities for the 157 Ta, 161 Re, 167 Ir and 185 Bi nuclei 
decaying from the 2s 1 / 2 state, the 53 9 Ise and 55 2 Css7 nuclei decaying from the lds/2 state, the gg 6 Tm77 
and 7 f 1 Lugo nuclei decaying from the O/in/2 state are predicted. Such spectra have orders of values similar 
Qjj ] to the experimental data for the bremsstrahlung photons emitted during the a-decay. This indicates on 

real possibility to study bremsstrahlung photons during proton decay experimentally and perform further 
d ' measurements. 

PACS numbers: 41.60.-m, 03.65.Xp, 23.50. +z, 23.20.Js 
Keywords: bremsstrahlung, proton-decay, tunneling, angular spectra 

q ■ 1 Introduction 

Last two decades many experimental and theoretical efforts have been made to investigate the nature of the 
bremsstrahlung emission accompanying a-decay of heavy nuclei [1T|22| and spontaneous fission |23H36j . A key 
idea of such researches consists in finding a new method of extraction of a new information about dynamics of the 
studied nuclear processes from measured bremsstrahlung spectra. Parallel study of the interested type of nuclear 
decay and the bremsstrahlung photons which are emitted during it, gives us more complete, richer information 
about the studied nuclear process. As examples, we established dependence between nuclear deformation and 
the bremsstrahlung spectrum accompanying a-decay of the 226 Ra nucleus [37j . Analysis of the bremsstrahlung 
emission which accompanies ternary fission of 252 Cf, allow us to understand better geometry of undergoing 
fission process in its first stage |38: ■ 

In our previous paper [22) we compared the bremsstrahlung spectra in a-decay of the 214 Po and 226 Ra nuclei, 
which have similar daughter nucleus-a particle potentials, and noted a clear difference between these spectra. 
In the cited paper we concluded that the different slopes of spectra were connected with the different Q- values of 
a-decay for two considered nuclei (Q a — 7.83 MeV for 214 Po and Q a = 4.87 MeV for 226 Ra) and we confirmed 
that such difference gives different contributions of emission from the tunnelling region into the total spectra. 
Now I put a question: which other parameters have important influence on the bremsstrahlung spectrum? 

Let us consider Figs. 4 and 7 in Ref. [32] where the calculated 7-ray emission probabilities for the spontaneous 
fission of 252 Cf are presented. One can see that the calculated spectrum is changed in dependence on mass split. 
So, one can suppose that the emission probability should be dependent on numbers of masses and charges of the 
daughter nucleus and the emitted fragment. The idea in Ref. [35] about influence of the electromagnetic charge 
of the daughter nucleus on the bremsstrahlung probability reflects this property only partially while effective 
charge of the decaying system is connected with such property more directly. Now if to consider the formula of 
the bremsstrahlung probability (for example, see (1) in |16j ) then one can find its direct dependence on square 
of the effective charge, Z^ s , i. e. we obtain a real basis for such supposition. 

As another type of decay which has larger effective charge in comparison with a-decay, the emission of 
proton from nucleus can be analyzed (for example, Z^ s equals to 0.286, 0.285, 0.288 and 0.303 for the proton 
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emitters 157 Ta, 161 Re, 167 Ir and 185 Bi, while it equals to 0.185 and 0.16 for the a-decaying nuclei 214 Po and 
210 Po, correspondingly). The first indication on possibility of presence of such process was given in [10] where 
authors estimated the averaged over all angles spectrum for the 113 Cs nucleus up to 250 keV of energy of the 
emitted photons. However, the detained study of such process has never been performed. First of all, a question 
is appearing whether the emission of photons accompanying such type of decay, is enough intensive to measure 
it experimentally (Q p -values for proton decay are smaller than for a-decay)? We are able to calculate the 
brcmsstrahlung probabilities for the a-decays, which turn out to be in enough well agreement with experimental 
data, without normalization of the theoretical spectra on experiment (see [1T1I42) ). Used in such calculations, the 
multipolc approach (started from [TTl[T5] ) looks to be the most accurate and corrected in angular description of 
the photons emission during the a-decay. It turns out that such model allows to calculate absolute values of the 
probability without any normalization relatively experimental data, and achieves enough good agreement with 
them. So, it could be a convenient basis for estimation of the probabilities of the bremsstrahlung photons in 
the proton-decay. According to our analysis, the dipole approximation of the wave function of photons (started 
from gives the overvalued spectra in absolute scale (see Fig. 9 in [13]), and by such a reason we shall not 
use it in this paper. 

Majority of nuclei emit protons in the state with nonzero orbital quantum number. So, from point of view 
of improvement of theory, it could be useful to generalize formalism of determination of matrix element of 
emission, taking into account such states. Inclusion of the states with nonzero orbital moment requires to add 
spin-orbital component into potential of interaction between the daughter nucleus and proton emitted (never 
considered before in this research). Note that up to-date it has not been known anything about the emission 
of photons from such states. So, it could be interesting to perform such calculations, that causes necessity of 
further improvement of theory of emission of bremsstrahlung photons in decays of nuclei. 

This paper answers on these questions, which is organized so. At first, the improved multipole model of the 
bremsstrahlung photons emitted during proton decay is presented, where emphasis is made on construction of 
the angular formalism of the matrix elements and calculation of the absolute bremsstrahlung probability. On 
its basis I perform theoretical study of the bremsstrahlung emission for some proton emitters, give predictions 
(in absolute scale) and analyze them. 

2 Model 

2.1 Matrix element of transition 

I define the matrix element like (2.11) in Ref. [15] (in the first correction of the non-stationary perturbation 
theory with stationary limits to = — oo and t\ = +00, and with normalization |C| — > 1): 



a fi = F fi ' 27T<5(U>/ - Wi + w), (1) 



where 



F fi = Zeff^yj^'pfakf), 



pfakf) = J2 e (a >'*p(fci,fc/), m 

1 9 \ £ ) 



p(ki,kf) 



and ipi(r) = \ki) and ipf( r ) = \kf) are stationary wave functions of the decaying system in the initial i-state 
and final /-state which do not contain number of photons emitted, Z e g and m are effective charge and reduced 
mass of this system. e( Q ) are unit vectors of polarization of the photon emitted, k is wave vector of the photon 
and w = k = |k|. Vectors are perpendicular to k in Coulomb calibration. We have two independent 
polarizations and e^ 2 ) for the photon with impulse k (a = 1,2). One can develop formalism simpler in the 
system of units where h = 1 and c = 1 , but we shall write constants Ti and c explicitly. Let us find also square 
of the matrix element a/, used in definition of probability of transition. Using the formula of power reduction 
of 8 -function (see [44], § 21, p. 169): 

[5(w)} 2 = 5{w) 6(0) = S(w) (2 7 r)~ 1 [ dt = S(w) (2^ T, (3) 
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we find (T — ¥ +00 is higher time limit): 

\a fl \ 2 =2ttT \F fl \ 2 ■ S(w f ~ w, + w), (4) 

that looks like (4.21) in Ref. [33] (with accuracy up to factor (2ir) 2 ) and like (42.5) in Ref. [JS] (exactly, see 
§ 42, p. 189). 

2.2 Linear and circular polarizations of the photon emitted 

Rewrite vectors of linear polarization e' a ' through vectors of circular polarization £ M with opposite directions 
of rotation (see Ref. [46], (2.39), p. 42): 

£-i = ^(e (1) -*e( 2 )), e+ i = -^(e«+ze( 2 )), & = e (3) = . (5) 
Then p(ki,kf) can be rewritten so: 

p(h,k f ) = Yl fr f (r)e- ikr ^-Mr)dv, (6) 

At=-i,i 

h± = T^, h- 1 + h +1 = -iy/2, £ eW^M-i+M+i. (7) 

V2 ce=l,2 

2.3 Expansion of the vector potential A by multipoles 

In further calculations of the matrix element p(ki,kf) the different expansions of function e~ ikr connected 
with the vector potential A of the electro-magnetic field of the daughter nucleus can be used. In spherically 
symmetric approximation I shall use the multipole expansion (see Ref. [46 , (2.106) in p. 58): 

^e lkr = M V27^(2/ + l) 1/2 i ' • [A ijU (r,M)+i M A^(r,£)], (8) 
l 

where (see Ref. (2.73) in p. 49, (2.80) in p. 51) 

A ijU (r,M) = ji(kr) T a ^(n ph ), 



21 + 1 Jl ~W > -"-H-1,M"P'W - ^ 



ji+i(kr) T u+1 ^(n ph ). 

Here, A/ A ,(r,M) and A/ M (r,i?) are magnetic and electric multipoles, ji(kr) are spherical Bessel functions of 
order I, T/;' iM (n) are vector spherical harmonics. We orientate the frame system so that axis z is parallel to the 
vector k. The functions T;// jM (n) have the following form (£0 = 0, see Ref. [35], p. 45): 

Tj/, m (n) = ^ (l,l,j\m- (i,n,m) Yi, m -^{n) (10) 
n=±i 

where (I, l,j | to— /z, n, m) are Clebsh-Gordon coefficients and Yi m (9, ip) are spherical functions defined, according 
to gS] (see p. 119, (28,7)-(28,8)). 

In the spherically symmetric approximation, wave functions of the decaying system in the initial and final 
states are separated into the radial and angular components, and these states are characterized by quantum 
numbers I and m. We shall be interesting in such photon emission when the system transits to superposition of 
all possible final states with different magnetic numbers m at the same orbital number I (for both states). So, 
let's write wave functions as 

ipi(r,li) = <pi(r,k) y^Yumj 

mi 

^/( r >*/) = <Pf(r,lf) 22 Y lfrn f (^{) 

m f 
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and obtain: 



where 



Pi i Pi 



Pi 



M 



and 



Pi 



ftp — 
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dr y dOr 2 V}(r) (j^w) A^foM), 
dr | dOr 2 ^(r) (|^(r)) Af„(r,.E). 



Using gradient formula (see (2.56), p. 46 in 



we obtain: 



where 



|U(r) = ^{^(r)5> iroj (<)} 



Zi / dtpi(r) k + l 



2k + l V dr 



k + l ( dtpi(r) k 



2k + 1 V dr 



(r)) ^T Mj+1 , mi (<), 



,1/ 



\J 2k + l Im ^ li ' l f' l P h ' li ~ 1 ) x 

^ 2Zf + 1 M ^ "'""'■) ' {"M^i ^/>'ph) — Zj ' "MZj, Z/, 



Zi (V + 1) 



V (2Z i + l)(2Zp h + l) MM/,U> 
x {jift,//,/ P h-l) + (Z* + 1)- J 2 (Z*, 

-^E (Zi, Z/, Zphj 



Zi ^ph 



V (2Zi + l)(2Z ph + l) 

x {Jl(ii, Z/, Zph + 1) + (k + 1) ' J2(k, 



+ 



I (k + m Ph + iT T n , . 

V(2/ i + l)(2/ ph + l) ' lE 

| Jl(Zj, If, Iph — 1) — Zi • hikjfjph — 1)| 



< - l,Z p h - l)x 
/^Ph-l)} - 
i - l,Z p h + l)x 

/^ph + i)} + 
i + l,Z p h - l)x 



(k + 1) z ph 



Ie (k,lf,lph,k + l,Z p h + l)x 



(2Z l + l)(2Z ph + l) 

x |Ji(Zi,Z/,Zph + i) - k ■ Mk,if,i P h + 1)}, 



Ji(k,lf,n) 
J 2 (k,lf,n) 



dr 



o 

+oo 



Im (h,lf,lph,h) = E EE* / ^("r^^^iWlTlip^Kt)^ 

Li™. m .. J 



rrii trif 



i E {h,i f MMM) = E EE^f^/K^y^^OTU^M^ 

Li™. ™ . */ 



(19) 



/Lf=d=l m/ 

For the case — we have essentially simpler formulas. In particular, for the matrix components we obtain: 



Pl P h - 


— Im (0, If, l t 


Pf = 


1 lph + 1 

V 2 u + 1 


+ 


/ l ph j 




V 2iph + 1 J 



(20) 



W o; P/ 'i i ^(0, if, iph, 1, lph + 1) ' J Qf,lpb + 1): 
W Zlph + 1 



where we use the simplified notation: 

J(l f ,n) = Ji(0,l f ,n). (21) 
Using the following value of the Clebsh-Gordon coefficient: 

(110 |1, -1, 0) = (110 |-1,1,0) = Jl, (22) 



from HUl and (HI) we obtain 



T 01 ,o(n;) = ( U0 I - Y h-M) ^ = Jl E £ 

f*=±l |U=±1 



(23) 



and for the angular integrals for transition into the superposition of all possible final /-states with different m / 
at the same I / from eq. (fT9"f we obtain: 

Iu(0,lf,l P h,n) = v/^ E E ( n > : Mpft|/ x_ /" t/ >M',A*) x 



lE{0,lf,l P h,l,n) = y ^ E ^ E ("^'^Im-m'jm'jM) x 
' JI=±1 Al'=±l 

x | ^(n^.-M'W^VM'^ 



(24) 



2.4 Vectors n* , n£, and calculations of the angular integrals at l { = 

Let us analyze physical sense of vectors n*, n£ and n p h- According to definition of wave functions ipi(r) and 
ipf(r), the vectors and n£ determine orientation of radius- vector r from the center of frame system to point 
where this wave functions describes the particle before and after the emission of photon. Such description 
of the particle has a probabilistic sense and is fulfilled over whole space. Change of direction of motion (or 
tunneling) of the particle in result of the photon emission can be characterized by change of quantum numbers 
I and m in the angular wave function: Yoo( n r) ~ * ^/m( n r) (which changes the probability of appearance of this 
particle along different directions, and angular asymmetry is appeared). The vector rip^ determines orientation 
of radius-vector r from the center of the frame system to point where wave function of photon describes its 
"appearance" . Using such a logic, we have: 

n ph = nl = nl = n r . (25) 

As we use the frame system where axis z is parallel to vector k of the photon emission, then dependent on r 
integrant function in the matrix element represents amplitude (its square is probability) of appearance of the 
particle at point r after emission of photon, if this photon has emitted along axis z. Then angle 9 (of vector 
n r ) is the angle between direction of the particle motion (with possible tunneling) and direction of the photon 
emission. 
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Let us consider the angular integral in (|24l) over dft. Using (j2"5")l . we find: 



3 (2^ + 1) (2n +1) (/ - 1)! (n-\m + //Q! x (2g) 



Y 32?r + (n+|m + /z'Q! 

y Pl (cos 6») (cos 6) Pl m+ v' I (cos 6>) • sin 6 dQ dip, 



where P; m (cos0) are associated Legandre's polynomial (see [45], p. 752-754, (c,l)-(c,4); also see [46] (2.6), p. 34) 
and the following restrictions on possible values of m and If have been obtained: 

m = -/z = ±l, Z/ > 1, n>\(x-(i'\ = \m + (x'\. (27) 

Let us introduce the following differential matrix elements dpf 1 and dpf dependent on the angle : 



dp M 



sin 6 d6 



I - ih+ l 



»^j(i f ,i Ph ) E mh -™ E (TsUj^CM 

m=±l fj,' = ±l 



d Pf = _ t h+i ph l l vh + l 



- J(lf,l ph - 1) 



sin 9 dO ]l 2l p h 

x E h ~™ E CyLu-i Ct-iW- (28) 



lf-\-l P h 



where 



E h -m E C l f d ph d ph + 1 flfd ph +l( )' 
m— ±1 ju'=±l 



C^Cln = (-l)^ +1 ^' + ^(n,l,« p ,|-m-/,',M',-m) 



'(2Z/ + 1) (2n+ 1) (n- |rn + /x'|)' 



32tt + (n+ |m + //])!' 

,|m+/i' 



(29) 



One can see that integration of functions (|28p by angle 6 1 with limits from to 7r gives the total matrix elements 
pf^ and exactly for transition into superposition of all possible final states with different m/ at the same If. 

We shall find the matrix element at the first values of If and l p h- We have If — 1, l p h = 1. Calculating 
coefficients and functions /-J^f* (6 1 ), from eq. (f2"5)) we obtain: 

^ - -| 4/1 • J(l,l)-sin 2 0cos0, 



sin 9 d6 8 V 7T 

sin 2 



= ^sin 2 .{V2.,(l,0) + ,(l,2).(l-3sin^)}. 



Integrating these expressions over angle 9, we find the integral matrix elements: 

Vf = 0, pf = iiyi.{j(l,0)-l^-J(l,2)}. (31) 

2.5 Angular probability of emission of photon with impulse k and polarization 

e («) 

I define the probability of transition of the system for time unit from the initial i-state into the final /-states, 
being in the given interval duf, with emission of photon with possible impulses inside the given interval dv p h, 
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so (see Ref. 05], (42,5) § 42, p. 189; Ref. [47 , § 44, p. 191): 



= -dv = 2-K \F Sl \ 2 6(wf -Wi+w)- dv, 



dv = dvf ■ dv p h, dv ph 



d 3 k w 2 dwd£l 



ph 



(32) 



(2tt) 3 (2ttc) e 



where dv v h and di// are intervals denned for photon and particle in the final /-state, d£l v h — d cos Q p h = 
sin 8 p h dOph dip p h, k p h — w/c. Ffi is integral over space with summation by quantum numbers of the system 
in the final /-state. Such procedure is averaging by these characteristics and, so, Ffi is independent on them. 
Interval dvf has only new characteristics and quantum numbers, by which integration and summation in Ffi 
was not performed. Integrating eq. (|32j) over dw and substituting eq. ([2]), we find: 



dW 



Z eff e2 hw fi 

m 2 2?r c 3 



p(h,kf) dVLp h dvf, Wfi = Wi — Wf = — l -— — -. (33) 



This is the probability of the photon emission with impulse k (and with averaging by polarization e( Q )) where 
the integration over angles of the particle motion after the photon emission has already fulfilled. 

I define the following probability of emission of photon with momentum when after such emission the 
particle moves (or tunnels) along direction n(: differential probability concerning angle is such a function, 
definite integral of which over the angle 9 with limits from to n equals to the total ■probability of the photon 



emission f33|). Let us consider function: 

dW(6 f ) z eff he2 w fi J n i ,dp*{ki,k f ,9 f ) 



{P(ki,kf) V" ,7 ' J; +h-e. ■ (34) 



dflphdcosOf 2tt c 3 m 2 \ ' d cos Of 

This probability is inversely proportional to normalized volume V . With a purpose to have the probability 
independent on V, I divide eq. (|34[) on flux j of outgoing a-particles, which is inversely proportional to this 
volume V also. Using quantum field theory approach (where v(p) = \p\/po at c = 1, see [44], § 21.4, p. 174): 

j = m v(pi), v, = |vi| = j?^ = - - ' kt , (35) 

tiii 

where rij is average number of particles in time unit before photon emission (we have = 1 for the normalized 
wave function in the initial i-state), v(pi) is module of velocity of outgoing particle in the frame system where 
colliding center is not moved, I obtain the differential absolute probability (while let's name dW as the relative 
probability): 

dP(<p f ,0f) dW(tp f ,6f) Ej 

dflph d cos Of dflph d cos Of hc 2 ki 

2 2 ( 36 ) 



2ir c 5 m 2 ki I d cos 



Note that alternative theoretical way for calculations of the angular bremsstrahlung probabilities in a-decays 
was developed in {20 based on different definition of the angular probability, different connection of the matrix 
element with the angle between fragment and photon emitted, application of some approximations. 

Let us find the bremsstrahlung probability at the first values U = 0, If = 1 and l p h = 1. Starting from 
eqs. (fT2|) and ([T3]) . and using the found differential and integral electrical and magnetic components (|30| and 
([3"T]) . I calculate: 



pi (h,kf) 



= -iJl'{j(l,0)-^V2-J(l,2)}, 

^ *f.{3J(l,l).cos*-V2J(l,0)- (37) 

sin dO 8 L 

- J(l, 2) ■ fl-3sin 2 0)\ -sin 2 
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and from eq. (|36[) I obtain the absolute angular probability: 



dP{ 



El+Ml 



(Of) 



d fl p h d cos 9 j 



Z eJf e Wfl Ei 



8nc 5 



.7(1,0) 
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V2- J(l,2) 



J*(1,0) + -^J*(l,2) • (l - 3 sin 2 0) 



V2 

4= J* (1,1) -cos ( 
V2 



(38) 



h.e. 



sin 



2.6 Proton— nucleus potential 



In order to describe the proton-decay, we need to determine the wave function of a nucleus consisting of a single 
proton interacting with a residual core. To describe the proton interacting with the spherically symmetric core 
we shall use the single-particle potential in standard optical model 



V(r) = v c (r) + v N {r) + v so (r) + Vi(r), 



(39) 



where v c {r), v^(r), v so {r) and vi(r) are Coulomb, nuclear, spin-orbital and centrifugal components given form 
in famous paper [35] (see eqs. (5) in the cited paper) as 



v N (r) 
v so (r) 
v c (r) = 



1 + exp 



r-R R 



vi{r) 



v; qi 

Ze^_ 

r 

2R r 



aR 
XI d_ 
dr 



1 + cxp 



1(1 + 1) 
2mr 2 

R 



t r - K so \ 



(40) 



3 - 



at r > R c , 
at r < R r - 



Here, Vr and V^ are strength of nuclear and spin-orbital components defined in MeV as (see eq. (8) in [4"8"] ) 

V R ; = 54.0 - 0.32 E + 0.4 Z/A 1 ? 3 + 24.0 7, V so = 6.2, (41) 

where 7 = (N — Z) /A, A and Z are mass and proton numbers of the daughter nucleus, E is incident lab energy. 
R c and 7?^ are Coulomb and nuclear radiuses of nucleus, aR and a so are diffusion parameters. In calculations 
I use (all parameters are defined in [48j ) : 



Rr^trA 1 / 3 , R c = r c A 1 / 3 , R^^r^A 1 / 3 , a R = 0.75 fm, 
tr = 1.17 fm, r c = 1.22 fm, r so = 1.01 fm, a so = 0.75 fm. 



(42) 



3 Results 

Let us estimate the bremsstrahlung probability accompanying the proton-decay. I calculate the bremsstrahhmg 
probability by eq. (J35J for k = or by eq. (l3l)|) for li ^ 0. The potential of interaction between the proton and 
the daughter nucleus is defined in eqs. (|39l) (|40|) with parameters calculated by eqs. (|4"Tj) - (j42"j) . To choose the 
convenient proton-emitters for calculations and analysis, I used systematics presented in Ref. (see Table 
II in the cited paper) and selected the 157 Ta, 161 Re, 167 Ir and 185 Bi nuclei decaying from the 2si/ 2 state (i.e. 
at k = 0), the 53 9 l56, 55 2 Cs57 nuclei decaying from the ld 5 / 2 state and the gg 6 Tni77, ^Lugo nuclei decaying 
from the 0/in/2 state (in order to analyze practical effectiveness of the model and convergence of calculations 
at li 7^ 0). The wave functions of the decaying system in the states before and after the photon emission are 
calculated concerning such potential in the spherically symmetric approximation. 

The result of calculations of the bremsstrahlung probabilities during proton decay of the nuclei pointed out 
above from from the initial state lsi/2 (at the chosen angle 9 = 90° between the directions of the proton motion 
(with its possible tunneling) and the photon emission) are presented in Fig. [TJ In Table [T] one can see some 
parameters for such studied nuclei. From here one can find that the different proton-emitters have practically 
similar effective charges, but different essentially tunneling regions. These difference between tunneling regions 
explain difference between the bremsstrahlung probabilities for such nuclei (that one can see in Fig.[TJ(a)). 



5 10 15 20 25 30 35 40 45 50 55 60 
40 80 120 160 200 240 280 320 360 400 4— ' — ' — ' — 1 — ' — <—•—> — ' — ' — ' — 1— — ' — ' — ' — ' — ' — ' — i — ■ 




40 80 120 160 200 240 280 320 360 400 5 10 15 20 25 30 35 40 45 50 55 60 



E., (keV) Radius, fm 

Figure 1: The bremsstrahlung during proton decay of the 157 Ta, 161 Re, 167 Ir and 185 Bi nuclei decaying from the 2si/2 
state (at the chosen angle 9 = 90°): (a) the absolute bremsstrahlung probabilities for these nuclei have the different 
trends that is explained by different tunneling regions for them, (b) the proton-nucleus potentials for these nuclei are 
similar practically. 



Proton decay data 


Parameters 


Turning points 


Tunneling 




Nucleus 


Q P , keV 


1 1/2 ' beC 


Rr, fm 


Vr, MeV 


1-st, fm 


2-nd, fm 


region, fm 


E ls 


73 ia 83 


947 


210 ms 


6.29 


-60.89 


7.25 


110.96 


103.71 


0.286259 


75 1 Re86 


1214 


180 


6.3517 


-60.8638 


7.32 


88.96 


81.64 


0.285328 


77 lr 90 


1086 


35 ms 


6.3912 


-61.0585 


7.32 


100.79 


93.46 


0.287924 


83 5Bi 98 


1611 


3.1 [AS 


6.6546 


-61.8599 


7.56 


74.28 


66.72 


0.303988 



Table 1: Parameters of the proton emitters decaying from the 2s!/ 2 state (E^ s is square of effective charge): 
one can see that these nuclei have similar effective charges but different Q-values and tunneling regions. This 
explains the different trends of the spectra presented in Fig. []] (a) . 
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Calculations of the bremsstrahlung probabilities for the case of Z, 7^ are essentially more complicated. 
Results of such calculations for the 5jj 9 l56, 55 2 Css7 nuclei decaying from the ld 5 / 2 state and the gg 6 Tm77, 
7f 1 Lu§o nuclei decaying from the Ohn/ 2 state are presented in next Fig. [5] (at angle 9 = 90°). One can see that 
these spectra have similar orders as results for the previous case at U — 0. According to analysis, inclusion of 




E(keV) EJkeV) 



Figure 2: (a) The bremsstrahlung accompanying proton decay of the 109 I and 112 Cs nuclei from the state ld$/2 (at 

8 = 90°); (b) the bremsstrahlung accompanying proton decay of the 146 Tm and 151 Lu nuclei from the state O/in/2 (at 

9 = 90°). 

the spin-orbital component v so into the total potential (|39p changes the resulting probability less than 1 percent 
(practically, 3-rd or 4-th digit in spectrum is changed; such spin-orbital influence on the spectrum has been 
estimated for the first time in the tasks of the bremsstrahlung during different types of nuclear decays). In next 
Fig. [3] angular distribution of the bremsstrahlung probability is presented. In particular, one can see that we 
obtain convergent calculations of angular integrals and the spectra at h =/= 0. 




Angle (degrees) Angle (degrees) 

Figure 3: The bremsstrahlung probability during proton decay of the 185 Bi (from the state 2si/ 2 ) and 109 I (from the 
state 1^5/2) nuclei in dependence on the 6 angle between direction of motion (with possible tunneling) of proton and 
direction of the emission of photon. In these figures different curves correspond to different energies of energy of the 
photon emitted. 



4 Conclusion 

In the paper the emission of the bremsstrahlung photons accompanying the proton decay of nuclei, is studied (for 
the first time, in systematic basis). The new improved multipolar model describing such a process is presented. 
The angular formalism of calculations of the matrix elements is stated in details. The absolute bremsstrahlung 
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probabilities for the 157 Ta, 161 Re, 167 Ir and 185 Bi nuclei decaying from the 2s 1 / 2 state (i.e. at U = 0), the 53 9 Is6 
and 55 2 Cs57 nuclei decaying from the 1^5/2 state (i.e. at k — 2), the gg 6 Tni77 and \\ Lugo nuclei decaying from 
the 0/in/2 state (i.e. at li = 5) are predicted. Such spectra have orders of values similar to the experimental 
spectra of the bremsstrahlung photons emitted during the a-decay (see [4][5l[T9l[2Tl[22] ) . It needs to note that 
arsenal of experimental measurements of the bremsstrahlung emission during different types of nuclear decays 
is not rich. A serious difficulty in obtaining of desirable accuracy in measurements lies in small values of the 
photon emission probability. One can hope that the presented theoretical estimations of the bremsstrahlung 
spectra indicate on the new possibility to perform further measurements and to study bremsstrahlung photons 
during proton decay experimentally. 



A Angular integrals at arbitrary ^ 

Let us calculate angular integrals: 

±1 nii rrif 

I E (k,lf,l P y,,h,h) = £ /^ ro /n/)T Wl , mi K)Tr ph i 2 ^(n ph )rfn, 

/a— ±1 nii rrif 

Substituting the functions Tj; im (n) in form (fTOj) . we obtain: 

Im (!i,lf,lph,h) = X] X! ^ X! {hA,h\mi- n',n',mi) x 

)i=ilmira; fj,'=±l 
X (Zph,l,Zph|M-M',M'>M) J Y l*m( n r) ■ ^i.nw-At'K) ' Y l* ph ^-^ i n ph) dfl, 

Ie (k,lf,l P h,h,h) = X! X! h > 1 X! {hA,h\mi- fi',fi',m t ) 

/i=±l mim/ fj, f = ±1 



(43) 



(44) 



x(l 2 ,l,l ph \fx-fx',fj,',fj,) J ^ m (nf)-yi lim( _^(n*)-y4 M _ M ,(n ph ) dfi. 

Here we have used the orthogonality condition for vectors £±1 (£0 = 0). 

Taking into account different variants of definition of spherical functions Yi rn (0, (p), we define them according 
to gS] (see. p. 119, (28,7)-(28,8)): 



■<»■*> " (- 1 »^ i '^™ P ' W( ™ 9) ' e '™' (45) 



Kir. 

where P™ (cos 0) are associated Legandre 's polynomials. Rewrite the angular integral in (|44|) as 



| Y* fnif (n r ) y/,., (n r ) y n * M _ M , (n r ) dfi = 



2l f + l(l f - \m f \)l hh + ljh-lmj- 2n+l(n-\n-n'\)\ > 

X V 4rr (Zy + lm/DlV 4tt (Z x + \rm - /*'!)! V 4 ^ + 1^ - ^'D' 



(46) 



i£7l /I 

J e i(_ m/+TOi _ M '_ M+ ^) v ^ . /" pl™/l(cos6») pl™'-^' 1 (cos 6>) P^-"'l(cosfl) • sin( 



2tt 

<y e H-tn/+"H-M- -*»+*«• ;v &p . J />'"'"(«.,* fc>) 7»'"-"'' Vosfl) /'./'•" lfrus«)-Hi,«,W. 


Integral over 95 in this expression is different from zero only in case 

fx = rrii — rtif. (47) 

So, we obtain the following restrictions: 

n > ll 1 ~ m'I = l m i ~ m / + m'Ii M — il- (48) 
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Using it, we find: 

/ Y i* m f( n r)Y ll:mi - tl '(n I )Y* tfl _ fl ,(n I ) dQ = 

— l_-[\lf+n+mi—n' ^i/+ii+n+|m/|+|m»— f/\+\rm— mj— fi'\ 



/ (2i / + l)(2fi + l)(2n+l) {l f -\m f \)l (h - |m. t - (n - |m, - m f - //|)! (49) 

X y 16tt (i/ + |m/|)! (Zi + \mi-n'\)l (n + |m<-m/ - fx'\)\ * 

7T 

x J P^ m/l (cos 0) pi™*""'! (cos 0) pl m '-™/-^l ( cos ) . sin 6 d6. 



Now we introduce the following generalized coefficient 



fjrrnnifn' _ f]\lf+l 2 +mi-vi' d f +h+h+\mf\+\mi-fi'\+\mi-mf-n'\ x 
hlflphhh v ' 

x (Zx, 1, I m» — //, //, mj) (Z 2 , 1, Z p h I TOj - m,f - //, fj,', rm - mf) x 



X 



' (2Zj + 1) (2fi + 1) (2Z 2 + 1) (// - K|)! (Zi-lm,-^!)! (fa - |m, - m f - »'\)\ 
16tt + (h + \mi-n'\)\ (h + \rm - m f - fi'\)l 



sin 6 d9 



sin 6» <Z0 

mimf fj,'=±l 
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and function 

y^m/M' (ff) = pl;»/l (cos0) J»l-HVl (00Bfl) pN*-™/-^l (00Bfl ). (51) 
In result, we obtain the following expressions for angular integrals Im and Ie- 

lM(li,lf,lpb,h)= E ( TOl -m f )h mi - mf E C hlZihl Ph ffl" W sin 0tZ0, 

rtiimf fi' — ±l q 

iE{k,i f ,i^,hM)= E E Cit / /W'w sin ^- 

m»ro/ ju'=±l q 

We also define differential functions form these angular integrals by angle so: 
d-Zjyf (Zi, Z/, Z p h, Zi) 
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